High Luminescence in Small Si/Si0 2 Nanocrystals: A Theoretical Study 
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In recent years many experiments have demonstrated the possibility to achieve efficient photo- 
luminescence from Si/SiC>2 nanocrystals. While it is widely known that only a minor portions of 
the nanocrystals in the samples contribute to the observed photoluminescence, the high complexity 
of the Si/Si02 interface and the dramatic sensitivity to the fabrication conditions make the iden- 
tification of the most active structures at the experimental level not a trivial task. Focusing on 
this aspect we have addressed the problem theoretically, by calculating the radiative recombination 
rates for different classes of Si-nanocrystals in the diameter range of 0.2-1.5 nm, in order to identify 
the best conditions for optical emission. We show that the recombination rates of hydrogenated 
nanocrystals follow the quantum confinement feature in which the nanocrystal diameter is the prin- 
cipal quantity in determining the system response. Interestingly, a completely different behavior 
emerges from the OH-terminated or Si02-embedded nanocrystals, where the number of oxygens at 
the interface seems intimately connected to the recombination rates, resulting the most important 
quantity for the characterization of the optical yield in such systems. Besides, additional conditions 
for the achievement of high rates are constituted by a high crystallinity of the nanocrystals and by 
high confinement energies (small diameters). 
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I. INTRODUCTION 



The indirect nature of the energy band-gap in silicon 
has always been the major obstacle for its employment 
in optical devices, since the momentum conservation 
requires additional mechanisms involved in the recom- 
bination process (e.g. electron-phonon interaction), 
that occur with low probability, hence producing very 
poor emitting rates. In the last decade the discovery of 
efficient visible photoluminescence (PL) and optical gain 
from silicon nanocrystals (NCs) have demonstrated the 
possibility to overcome the indirect band gap of silicon 
by exploiting the quantistic behaviour of the matter at 
the nanoscale. 1,2 Theoretically, the optical emission has 
been attributed to transitions between states localized 
inside the nanocrystal [as a consequence of the so-called 
quantum confinement (QC) effect], 3-6 or between defect 
states. 7-13 While there is still some debate on which of 
the above mechanisms primarily determines the emission 
energy, some works have proposed that a concomitance 
of both mechanisms is always present, favouring one or 
the other depending on the structural conditions. 14-20 
Embedding Si-NCs in wide band-gap insulators is 
one way to obtain a strong QC. Si-NCs embedded 
in a silica matrix have been obtained with several 
techniques as ion implantation, 9,21 chemical vapour 
deposition, 13 ' 20,22 ' 23 laser pyrolysis, 8,24 electron beam 
lithography, 25 sputtering, 10,26 and some others. Inde- 
pendently on the fabrication technique, in experimental 
samples no two NCs are the same, practically binding 
the measures to collective properties. Experimentally, 
several factors contribute to make the interpretation 
of measurements on these systems a difficult task. 



For instance, samples show a strong dispersion in the 
NC size, that is difficult to be determined. In this 
case it is possible that the observed quantity does not 
correspond exactly to the mean size but instead to 
the most responsive NCs. 27 Again, NCs synthesized by 
using different techniques often show different properties 
in size, shape and in the interface structure. Also, in 
solid nanocrystal arrays some collective effects caused 
by electron, photon, and phonon transfer between the 
dots can strongly influence the luminescence dynamics 
in comparison with the case of isolated NCs. 9 
Recently, some studies on the PL spectra of individual 
silicon quantum dots have reported lincwidths of 2 
meV at T~35 K, clearly below k B T at this temper- 
ature, demonstrating true quantum dot PL emission 
characteristics. 24 ' 25,28 The treatment of single NCs 
with atomic-like models allows a comparison between 
the collective response of a sample with the sum of 
the responses of the individual NCs, leading to useful 
informations about the NC-NC interaction mecha- 
nisms. Also, many works have already indicated the 
predominant contribution from Si-NC with particular 
size/shape, 27,29,30 providing important implications for 
the identification of the best conditions for light emission 
and optical gain. 



In this work we calculate the recombination-rates (RRs) 
of several Si-NCs in order to provide a description of the 
optical emission yield of such systems. Then, a system- 
atic investigation on the different conditions of size, oxi- 
dation, strain, and amorphization is performed, allowing 
the identification of the circumstances that maximize the 
luminescence. In Sec. II the structures and the methods 
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TABLE I: Structural characteristics for the crystalline embedded NCs (top set), amorphous embedded NCs (center set), crys- 
talline suspended NCs (bottom set). For each structure are reported, respectively: number of NC silicons, number of core silicons 
(not bonded with oxygens), centered or not on one silicon, number of oxygens bonded to the NC, number of oxygens bridging 
two NC silicons, average diameter d, supercell volume V s . 



used to calculate the emission rates are described. The 
results are presented and discussed in Sec. III. Conclu- 
sions are reported in Sec. IV. 

II. STRUCTURES AND METHOD 

We built three different sets of NCs (sec Table I): 
crystalline NCs embedded in Si02 (generated from 
betacristobalitc), amorphous NCs embedded in silica 
(generated from a glass), OH-terminatcd NCs suspended 
in vacuum. A further set of hydrogenated NCs is 
obtained from the latter by replacing the OH groups 
with hydrogens and a subsequent geometric relaxation. 
The crystalline embedded structures have been obtained 
from a betacristobalitc cubic matrix by removing all 
the oxygens included in a cutoff-sphere, whose radius 
determines the size of the NC. By centering the cutoff- 
sphere on one silicon or in an interstitial position it is 
possible to obtain structures with different symmetries. 
To guarantee a proper shielding of the introduced strain 
by the surrounding silica we preserved a separation of 
about 1 nm between the NCs replica. 
The glass models have been generated using classical 
molecular dynamics simulations of quenching from a 
melt, followed by ab-initio relaxations (See Ref. 31 for 
further details). The amorphous dot structures have 
been obtained, as for the crystalline systems, by applying 
the cutoff sphere on the glass supcrcclls. It is worth 
to note that also in the crystalline case the embedding 
matrix gets amorphized after the inclusion of the NC, 



due to the metastable nature of the betacristobalite. 31 
The hydroxided NCs have been obtained by cutting 
off a spherical region from a bulk-silicon supercell, and 
successively passivating the dangling bonds with OH 
groups. 

The relaxation of all the structures have been performed 
using the SIESTA code 32,33 and Troullier-Martins pseu- 
dopotcntials with non-linear core corrections. A cutoff 
of 150 Ry on the density and no additional external 
pressure or stress were applied. Atomic positions and 
cell parameters have been left totally free to move. 

Each set is formed by NCs of different size, inter- 
face configuration, and symmetry. While the variation 
of the size is straightforwardly connected to the con- 
finement energy, the latter aspects are associated to the 
formation of defect-states and to the localization of the 
electron- hole pairs. 

For each structure we calculated the eigenvalues and 
cigenfunctions of the independent-particle Hamiltonian, 
using the Density Functional Theory (DFT), with 
the ESPRESSO package. 34 Calculations have been 
performed using norm-conserving pscudopotcntials 
within the Local Density Approximation (LDA) with 
a Cepcrlcy- Alder exchange-correlation potential, as 
parametrized by Perdew-Zunger. An energy cutoff of 60 
Ry on the plane wave basis have been considered. 
It is well known that the DFT-LDA severely underesti- 
mates the band gaps for semiconductors and insulators: 
since DFT is a ground state theory, the electronic 
energy eigenvalues do not correctly describe the ex- 
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citation energies of electron and holes. A correction 
to the fundamental band gap is usually obtained by 
calculating the separate electron and hole quasiparticle 
energies via the GW method. 35 The true quasiparticle 
energies, however, are still not sufficient to correctly 
describe a process in which electron-hole pairs are 
created. Their interaction can lead to a dramatic shift 
of peak positions as well as to distortions of the spectral 
lincshapc. Besides, previous many-body calculations on 
Si-NCs reported absorption spectra very close to the 
DFT-RPA calculated ones. 31 ' 36 ' 37 In Rcf. 31 we verified 
this statement for the Siio and a-Siio embedded NCs, 
showing that the self-energy (calculated through the 
GW method) and electron-hole Coulomb (calculated 
through the Bcthc-Salpeter equation) corrections exactly 
cancel out each other (with a total correction to the 
gap smaller than 0.1 cV) when the local field effects 
(LFE) are neglected. Furthermore, some our recent 
calculations (still unpublished), show that the LFE 
actually blue-shifts the absorption spectrum of the 
smallest systems (d<l nm), with corrections of the order 
of few tenths of eV. Instead, for larger NCs no blue-shift 
is observed. Therefore, while such correction should be 
taken into account for a rigorous calculation of the RRs, 
we expect that the LFE would affect the final rates by 
only a few percent. Moreover, an enlargement of the gap 
for the smallest NCs would increase the calculated RRs, 
strengthening our conclusions about the size effects. 
The reasons discussed above justify the choice of DFT- 
LDA for our calculations, allowing a good compromise 
between results accuracy and computational effort. 

For a first-order radiative process, the spontaneous decay 
probability for unit time from a state \i) to a state \j), 
with Eij = Ei — Ej > 0, is defined from the Fermi golden 
rule and is given by 38 

l6TT 2 ne 2 Eij 2 
3 m W MP'*)' ' (1) 

where n,e,m,h,c,p are the refractive index, electron 
charge, electron mass, Planck constant, light speed, and 
momentum operator, respectively. The static refractive 
index has been obtained from the real part of the di- 
electric function at zero energy, n = e\[oj = 0), cal- 
culated in the random-phase approximation (RPA). We 
observed an increasing trend of n with the NC diameter, 
with 1.1 < n < 1.2 (1.75 < n < 1.85) for the suspended 
(embedded) NC sets. Considering the large variation of 
the final emission rates (practically determined by the 
matrix elements) we approximated n = 1 for all the cal- 
culations. Equation 1 does not take into account any 
summation over the k-points of the reciprocal space, since 
the states belonging to the NC are strogly localized (due 
to the high QC), and no energy dispersion in the bands 
occurs (see Ref. 31). For this reason we performed the 
calculation of the probabilites at the single V point. 
From Eq. 1 we can define the thermally averaged recom- 



bination rate by 39 : 

where i and j are states belonging to the conduction and 
valence bands, respectively. The spin is not included in 
the RR calculation, and therefore the effects due to ex- 
change interaction arc neglected, with a possible influ- 
ence on the low temperature (< 50 K) dependence of the 
calculated rates. 40 Morerover, while the recombination 
through nonradiative decay mechanisms is expected to be 
favoured for high temperatures, all experimental results 
indicate that radiative transitions dominate for a broad 
temperature range, up to room temperature. 41 ' 42 Finally, 
Eq. 2 assumes that the thermalization of the electron and 
the hole after the excitation in the bands is more efficient 
than the radiative recombination. This condition is rea- 
sonably satisfied because the radiative recombination is 
usually not very efficient: RRs for spherical crystallites 
of different diameters have been reported ranging from 
about to 0.1 ns _1 , 39 compared to timescales of the 
order of 1 ps for thermally-activated processes. 

III. RESULTS 

In this Section we present the thermally averaged 
RRs, calculated through Eq. 2, for the whole set of 
NCs. Figure 1 reports the RRs at the temperatures 
of 5 K, 300 K, and 600 K, as a function of the energy 
gap, Eq [given by the difference between the LUMO 
(Lowest Unoccupied Molecular Orbital) and the HOMO 
(Highest Occupied Molecular Orbital) states], and of 
the oxidation (or hydrogenation) degree, Q [given by 
the ratio between the oxygens (or hydrogens) at the 
interface by the number of silicons to which they are 
connected], for the OH-terminated and H-terminated 
NCs. The Eq and f2 values for each suspended NC 
are reported in Table II. Some trends clearly emerge 
from the rates reported in Fig. 1. The simplest case of 
hydrogenated NCs shows a clear increase of the emission 
rate with the gap (Fig. la), indicating that the smaller 
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TABLE II: Oxydation degree ft, and HOMO-LUMO gaps for 
the H- and OH-terminated NCs. 
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FIG. 1: Thermally averaged recombination rates, calculated 
at different temperatures, for H-terminated [a) and b)J and 
OH-terminated [c) and d)] NCs, as a function of the HOMO- 
L UMO energy gap or of the oxidation /hydrogenation degree, 
fi. The dashed line stands for the power-law fit to the data at 
T = 300 K (see text). 



NCs are the more efficient, with characteristic times 
even below 0.1 /is. A power-law fit to the data at 



T = 300 K results 7.7 • 10 6 (£G/eV-2.47) 



1.19 



(Fig. 



la, dashed line). By considering the proportionality 
introduced by the Eij term in Eq. 1, the fit points out 
an almost constant matrix element with the NC size, 
as occurring for an ideal two-level exciton system. 42 
Besides the general trend, some scattering is present in 
the data, that can be better interpreted by examining 
the RR-vs-fi dependence (Fig. lb). It is well known 
that the HOMO-LUMO gap of hydrogenated NCs is 
not related to the passivation regime, simply following 
the QC. 43 This uncorrelation is well reflected in the 
RR-vs-fi plot, showing a much clearer picture in that 
most of the NCs emits at the same rate (of about 10' 
s _1 ) independently of 51. Anyway, for NCs with the 
same f2 the characteristic time increases with size, as 
shown for f2~1.5 and 51=3. We note that the scattering 
of the data is more pronounced at low temperatures, 
where the RR is mostly due to the LUMO-HOMO 
recombination. Such behavior is due to the fact that 
the lowest transition corresponds rarely to the most 
efficient among all the possible transitions. Contrary, 
at high temperatures more channels contribute to the 
average, leading to a smoother trend of the RR. Quan- 
titatively, the variation of the RR in the temperature 
range of 300-600 K is always smaller than one order 
of magnitude, making such systems poorly dependent 
on T. The above results for H-tcrminated NCs are 
in good agreement with calculations on porous silicon 
in the tight-binding approximation, 39 and with other 
DFT-based calculations on Si/Gc NCs. 44 
The analysis of the data relative to the hydroxided 
NCs follows that of the hydrogenated ones. The RR 
dependence on Eg (Fig. lc) shows a supra-linear trend, 
with some noticeable deviations induced by the Ec-fl 
correlation. The power-law fit to the data at T = 300 
K results 3.7 • 10 5 (.E G /eV-0.71) 2 - 94 s"\ very near 
to the cubic relation occurring for an ideal two-level 
atom. 42 Besides the QC effect, clearly dominating for 
hydrogenated NCs, in the case of hydroxided NCs we 
have to consider that the band-gap is strongly correlated 
to the oxidation degree, producing strong fluctuations 
of the gap with the NC size. 43 Such correlation could 
be responsible for the large broadening of the PL 
spectra that has been observed experimentally, also at 
low temperatures, 3,7, 10,19-24,45-47 even for apparently 
monodispersed multilayered samples. 49 This consider- 
ation reflects on the RR-vs-£l dependence (Fig. Id). 
Contrary to the hydrogenated case, here some depen- 
dency on fi of the RR occurs, with some correlation 
emerging between the two quantities. In summary, we 
observe that the hydrogenated NCs present higher rates 
with respect to OH-terminated NCs, as experimentally 
demonstrated. 47,50 Besides, for the hydroxided NCs the 
optical emission is favoured when a high O/Si ratio is 
achieved, in well agreement with some experimental 
observations. 16,48 
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The calculation of the RRs for the embedded structures, 
in comparison with those of the hydroxidcd-suspendcd 
NCs, could shed some light on the role of the embedding 
matrix. For example, it is known that the strain induced 
by the silica on the NC tends to red-shift the band-gap 
of the smallest NCs (< 1 nm) while blue-shifting that of 
the largest ones (> 2 nm). 43,51 ' 52 Therefore, the effects 
of such strain on the emission properties are investigated 
by comparing the results for the suspended (unstrained) 
and the embedded systems. Figure 2 reports the RRs for 
the embedded NCs set, at the temperatures of 5 K, 300 
K, and 600 K, as a function of the HOMO-LUMO gap 
(top panel), and of fi (bottom panel). In this case ft is 
also an index for the strain induced on the NC by the 
silica matrix, as discussed in Rcf. 43 where a correlation 
between the two has been proposed. The calculations for 
the embedded structures confirm an increase of the rates 
with E G and 0. The fit to the data at T = 300 K re- 
sults 4.0 • 10 5 (£ G /eV-0.45) 2 - 94 s -1 , very close to what 
obtained for the hydroxided NCs. Therefore, the strain 
seems to not modify the picture, discussed above, regard- 
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FIG. 2: Thermally averaged recombination rates, calculated 
at different temperatures, for the Si/SiOi embedded NCs, a) 
as a function of the HOMO-LUMO energy gap, and b) as a 
function of the oxidation degree, Q,. The dashed line stands 
for the power-law fit to the data at T = 300 K (see text). 



ing the effects of size and oxidation. The two smallest, 
highly oxidized NCs (S12 and Sis) follows the trend, at- 
taining impressive characteristic times, even below 1 ns 
(Sis case). Nevertheless, the embedding matrix has some 
effects on the final rates, that are sometimes strongly 
increased (Siu), sometimes unaffected (Si32), and other 
times strongly reduced (Siio and Si35). From the compar- 
ison of Fig. Id and Fig. 2b, we note that the introduction 
of the matrix produces a clear RR-vs-Sl trend, revealing 
that the oxidation is a crucial quantity for the character- 
ization of the optical properties of Si/Si02 systems. 
The opportunity of considering two opposite ideal cases 
of cristalline and amorphous NCs comes from the fact 
that real samples are characterized by a certain amount 
of amorphization, especially for smaller NCs. 17,53 ' 54 Im- 
portant connections of the amorphization with the PL 
spectrum have been reported, generally indicating, for 
amorphous NCs, reduced 55,56 or comparable 45 ' 53,57 PL 
intensities with respect to the cristalline case. 
Our results show that the introduction of disorder seems 
to flatten the RR trend, with all the values oscillating 
within 10 6 -10 7 s" 1 . We note that the higher emission 
of the a-Siio NC with respect to the Siio case can be 
explained with the larger oxidation degree of the amor- 
phized NCs (see bottom panel of Fig. 2). The same rule 
applies to the a-Si32 NC, presenting a smaller RR due to 
a reduced oxidation with respect to its crystalline coun- 
terpart. Instead, when the oxidation degree is preserved, 
the disorder always reduces the RR (Si2, Sis). Therefore, 
our results confirms the experimental observations, indi- 
cating that the amorphization process generally reduces 
the optical yield. 



IV. CONCLUSIONS 

We evaluated the recombination rates at different tem- 
peratures for Si-NCs of diameters in the range of 0.2-1.5 
nm, as for H- and OH-terminated NCs, like as for NCs 
embedded in a Si02 matrix. Different conditions of size, 
passivation, strain, and symmetry have been considered 
in order to investigate for the best conditions of radiative 
emission. 

We found that the rates follow a trend with the emission 
energy that is nearly linear for the hydrogenated NCs, 
and nearly cubic for the NCs passivated with OH groups 
or embedded in Si02- Moreover, the hydrogenic passi- 
vation produces higher optical yields with respect to the 
hydroxilic one, as evidenced experimentally. 47,50 Besides, 
for the hydroxided NCs the emission is favoured for sys- 
tems with a high O/Si ratio. 16 ' 48 

The analysis of the results for the embedded NCs re- 
veals a clear picture in which the smallest, highly oxi- 
dized, crystalline NCs, belong to the class of the most 
optically-active Si/Si02 structures, attaining impressive 
rates of less than 1 ns, in nice agreement with experimen- 
tal observations. 45 ' 53, 55-57 
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